The Melting Curve and Premelting of MgO

The melting curve for MgO was obtained using molecular dynamics and a non-empirical, many-body potential. We also studied premelting effects by computing the dynamical structure factor in the crystal on approach to melting. The melting curve simulations were performed with periodic boundary conditions with cells up to 512 atoms using the ab-initio Variational Induced Breathing (VIB) model. The melting curve was obtained by computing ∆H m and ∆V m and integrating the Clapeyron equation. Our ∆H m is in agreement with previous estimates and we obtain a reasonable ∆V m , but our melting slope dT/dP (114 K/GPa) is three times greater than that of Zerr and Boehler [1994] (35 K/GPa), suggesting a problem with the experimental melting curve, or an indication of exotic, non-ionic behavior of MgO liquid. We computed S(q, ω) from simulations of 1000 atom clusters using the Potential Induced Breathing (PIB) model. A low frequency peak in the dynamical structure factor S(q, ω) arises below the melting point which appears to be related to the onset of bulk many-atom diffusive exchanges. These exchanges may help destabilize the crystalline state and be related to intrinsic crystalline instability suggested in earlier simulations.